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ABSTRACT 

This  paper  focuses  on  potential  false  conclusions  regarding  polarization  parameters  obtained  through  the  use  of  a  digital 
camera.  It  outlines  a  methodology  for  calibrating,  displaying  and  visualizing  polarization  profiles  of  natural,  daylight  scenes 
using  a  COTS  (commercial-off-the-shelf)  digital  camera.  It  presents  cause  and  effect  scenarios  related  to  spatial  and 
temporal  misregistration  of  polarization  data  and  the  acquisition  of  accurate  polarization  parameters,  such  as  the  degree  of 
polarization  and  the  polarization  azimuth  and  ellipticity  angles.  Trade-offs  between  precise  spatial  registration  and  precise 
temporal  registration  are  discussed. 

THE  STOKES  PARAMETERS 
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The  human  visual  system  can  barely  detect  polarized  light.  Some 
observers  can  determine,  through  sighting  Haidinger’s  brush1,  that 
skylight  is  polarized.  However,  the  human  visual  system  is  incapable 
of  determining  the  complete  polarization  state  of  a  light  beam.  To 
determine  the  state  of  polarization  corresponding  to  points  in  a  scene, 
three  independent  parameters  must  be  determined  for  each  point2. 

For  example,  the  three  independent  parameters  could  be  file 
amplitudes  of  the  x  and  y  components  of  the  electric  vector  Eox  and 
Eoy  and  their  phase  difference,  8,  along  the  optical  axis.  A  prominent 
method  to  determine  these  three  independent  parameters,  using 
measurable  quantities,  is  the  Stokes  method3.  This  method  involves 
measuring  four  intensities  of  a  light  beam.  Each  measurement 
corresponds  to  the  intensity  of  the  beam  after  it  passes  through  each 
of  four  different  filter  system  arrangements.  The  four  Stokes 
parameters,  sometimes  called  So,  Si,  S2  and  S3,  are  derived  from 
these  four  measured  intensities  and  form  a  four-element  column 
vector  in  four-dimensional  mathematical  space.  The  Stokes 
parameters  are  applicable  to  a  beam  of  light -that  is  completely 
polarized,  partially  polarized  or  unpolarized;  the  beam  may  be 
monochromatic  or  polychromatic.  The  Stokes  parameters  for  completely  polarized  light  propagating  along  the  +z-axis  are 


Figure  1  The  relationship  between  the  Stokes 
parameters  and  the  Poncard  sphere. 


So 13  Eox2  +  Boy2  Si  =  Eox2  —  Eoy2  Sj  =  2  Eox  Eoy  cos  8  S3 =  2  Eox  Eoy  sin  8 


(1) 


Normalized  parameters  are  obtained  by  dividing  So,  Si,  S2  and  S3  by  So. 
polarization  P  is  given  by4 


For  partially  polarized  light  the  degree  of 
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where  Ipotn-tn*)  is  the  intensity  of  die  polarized  component  and  ITolai  is 
the  total  intensity.  When  the  Stokes  parameters  are  normalized,  the 
degree  of  polarization  P  becomes  the  radius  of  a  sphere  with  center 
coincident  to  the  typical  Poincafo  sphere  (see  Fig.  1)  of  unit  radius  (P 
=  1).  The  inside  of  die  Poincard  sphere  can  be  used  to  represent  the 
polarized  portion  of  partially  polarized  light  (P  <  1).  The 
polarization  forms  on  the  surface  of  the  inside  spheres  (P  <  1)  are 
exactly  the  same  as  the  polarization  forms  on  the  surface  of  the 
sphere  of  unit  radius.  From  the  geometry  of  Figure  1  and  Eq.  (2)  it 
can  be  shown  that  for  normalized  Stokes  parameters  and  for  any 
degree  of  polarization 
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Si  =  P  cos  2\p  cos  2%  S2  =  P  sin  2vp  cos  2%  S3  =  P  sin  2x  (3) 

Figure  2  The  polarization  ellipse.  The 

In  fact,  Eq.  (3)  represents  the  spherical  coordinates  for  any  point  on  parameters  y,  x  are  the  polarization  azimuth  and 
the  surface  or  inside  of  die  Poincard  sphere  where  (x,  y,  z)  =  (Si,  S2,  ellipticity  angles. 

S3). 


The  polarization  form  associated  with  elliptical  polarization  is  given  by  the  following  expression,  as  shown  in  Fig.  2 
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The  polarization  azimuth  angle,  ip,  and  die  ellipticity  angle,  x,  are  also  defined  in  Fig.  2.  Using  Eq.  (3),  ip  and  %  can  be 
determined  from 

sin  2z  =  -j ■  and  tan2*/=|-2-  (5) 

VS^+S^+Sj2  s, 

The  parameter  x  varies  from  +45°  to  -45°;  it  is  positive  for  right-handed  polarization  forms  and  negative  for  left-handed 
polarization  forms.  The  parameter  ip  varies  from  0°  to  180®;  it  is  0°  for  horizontal  polarization  forms  and  90°  for  vertical 
polarization  forms. 

MEASURING  THE  STOKES  PARAMETERS 


The  Stokes  parameters  can  be  measured  using  a  linear  polarizer  and  a  retarder.  If  light  is  first  transmitted  through  a  linear 
retarder  with  its  fast  axis  oriented  at  an  angle  Q  with  respect  to  the  x-axis  and  then  transmitted  through  a  linear  polarizer  with 
its  transmission  axis  oriented  at  an  angle  0  with  respect  to  the  x-axis,  die  intensity  of  the  transmitted  light  can  be  expressed  as 

I(Q,e,e)  = 

Eox2  [  cos2  0  -  sin  n  cos  n  sin  2(Q  -  0)  (1  -  cos  e  )]  +  Eoy2  [sin2  0  +  sin  Q  cos  O  sin  2(Q  -  0)  (1  -  cos  e  )]  + 

Eox  Eoy  [  cos  2(Q-0)  sin  2Q  cos  8  -  sin  2(Q-0)  (cos  2D  cos  e  cos  8  +  sin  s  sin  8)  ]  (6) 

where  8  *  (<(ix  -  4>y)  is  the  phase  difference  between  the  x  and  y  components  of  the  incident  light  and  e  is  the  phase  difference 
produced  by  the  retarder.  I(Q,0,e)  denotes  an  intensity  measurement  corresponding  to  a  particular  set  of  values  for  Q,  0  and 
e.  For  a  quartz  retarder5 


-50.876^  (7) 
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where  XT  is  the  wavelength  for  6  =  90°  (the  tuned  wavelength)  and  A,  is  the  incident  wavelength,  both  in  nanometers.  Figure  3 
shows  the  phase  difference  for  a  quartz  retarder  as  a  function  of  the  input  wavelength  for  selected  tuned  wavelengths.  The 
smallest  variation  in  e  occurs  for  a  quartz  retarder  tuned  to  the  blue  end  of  the  visible  spectrum  and  for  input  wavelengths  in 
die  red  end  of  the  visible  spectrum  (die  red  channel  of  a  COTS  digital  camera). 


The  Stokes  parameters  can  be  determined  experimentally  from 
the  following  four  measurements,  as  shown  in  Fig.  4: 

I,  =  1(0, 0,0)  I2  =  1(0, 90,0) 

I3  =  I(0,45,0)  14  =  1(0,45,8)  (8) 

In  this  arrangement,  a  retarder  is  inserted  into  the  optical  path  for 
only  the  fourth  measurement.  The  four  measurements  described 
in  Eq.  (8)  are  not  the  only  four  measurements  that  will  lead  to  a 
determination  of  the  Stokes  parameter.  In  a  previous  work5,  we 
describe  a  simple  two-element  variable  filter  system  where 
neither  the  polarizer  nor  the  retarder  needs  to  be  removed  for  any 
measurement. 

Figure  5  shows  an  experimental  setup  to  acquire  this  four-image 
sequence.  It  shows  a  digital  camera,  a  linear  polarizer,  a  hinged 
retarder  and  a  hinged  filter  mounted  on  a  tripod.  The  retarder 
and  filter  are  hinged  so  that  they  can  be  moved  quickly  and 
precisely  into  and  out  of  the  optical  path  perpendicular  to  the 
optical  axis.  Substituting  Eq.  (8)  into  Eq.  (6)  and  comparing  the 
results  with  Eq.  (1)  confers  that 

So  =  1(0, 0, 0)  + 1(0, 90, 0)  =  I,  + 12 

S,  =  1(0, 0, 0)  - 1(0, 90, 0)  =  I,  - 12 ; 

S2  =  2 1(0, 45, 0)  -  S0=  2 13  -  So 


S3  = 
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sine 


(9) 


Quartz  Retarder  Phase  Difference 


Input  Wavelength  (nm) 


Figure  3.  Phase  difference  produced  by  a  quartz 
retarder  tuned  to  selected  wavelengths. 


Il  =  1(0,0, 0)  12  =  1(0,90,0)  fa  - 1(0,45,0)  U  =  1(0,45,0) 


Figure  4  The  four-filter  system  used  to  determine  the  Stokes 
parameters.  P  is  the  transmission  axis  of  a  linear  polarizer 
and  R  is  the  fast  axis  of  a  zero-order  quartz  quarter- wave 
retarder. 

Since  the  Stokes  parameters  require  quantities  proportional  to 
light  intensity,  it  is  essential  that  the  digital  camera  be  calibrated 
to  convert  RGB  pixel  values  to  intensities.  This  task  is  easily 
accomplished  by  photographing  a  white  standard  illuminated  by 
a  collimated  white  beam  of  light  with  a  large  number  of  different 
neutral  density  filters  in  front  of  the  digital  camera.  The  RGB  to 
optical  density  conversion  equation  for  an  Epson  850Z  digital 
camera  was  determined  to  be  of  the  form 


Figures  An  experimental  camera  setup  to  acquire 
the  four-image  sequence  for  obtaining  the  Stokes 
parameters  associated  with  each  pixel  in  a  scene. 


Ax2  +  Bxy  +  Cy2  +  Dx  +  Ey  +  F  =  0  (10) 

where  x  and  y  correspond  to  optical  densities  and  RGB  values 


respectfully.  The  coefficients  in  Eq.  (10)  are  different  for  each  color  channel  of  the  digital  camera.  Using  the  fixed  white 
balance  option  of  the  Epson  850Z,  the  coefficients  obtained  for  the  green  channel  are: 
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This  paper  has  outlined  a  methodology  for  measuring  the  four 
Stokes  parameters  and  calibrating,  displaying  and  visualizing 
polarization  profiles  of  natural,  daylight  scenes  using  a  COTS 
digital  camera.  Measurement  errors  due  to  spatial  and  temporal 
misregistration  of  polarization  data  have  been  described  and 
examples  of  each  have  been  given.  It  has  been  shown  that  when 
the  intensity  increases  between  image-acquisition,  the  x-values 
increase  and  produce  pseudo  right-handed  elliptically  polarized 
light.  When  the  intensity  decreases  between  image-acquisition, 
the  x-values  decrease  and  produce  pseudo  left-handed  elliptically 
polarized  light.  Trade-offs  between  precise  spatial  registration 
and  precise  temporal  registration  have  also  been  discussed. 
Current  studies  strive  to  reduce  the  time  required  to  obtain  the 
four  Stokes  parameters  and  to  minimize  spatial  and  temporal 
misregistration. 


Figure  15.  An  apparatus  that  simultaneously  creates 
four  different  attenuated  images  of  a  scene. 
Corresponding  pixel-values  from  each  image  are 
used  to  calculate  the  Stokes  parameters  for  each  pixel 
location. 
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